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Dark matter within high surface brightness spiral galaxies 

Thilo Kranz^'^, Adrianne Slyz^ and Hans- Walter Rix^ 

§ : ABSTRACT 

(N : 

Q I We present results from a detailed dynamical analysis of five high surface 

Q ; brightness, late type spiral galaxies NGC 3810, NGC 3893, NGC 4254, NGC 

CN I 5676, and NGC 6643 which were studied with the aim to quantify the luminous- 

I to-dark matter ratio inside their optical radii. The galaxies' stellar light distri- 

^ I bution and gas kinematics have been observed and compared to hydrodynamic 

O I gas simulations which predict the gas dynamics arising in response to empiri- 

^ I cal gravitational potentials, which are combinations of differing stellar disk and 

CN ' dark halo contributions. The gravitational potential of the stellar disk was de- 

^ ' rived from near-infrared photometry, color-corrected to yield a constant stellar 

O ■ mass-to-light ratio (M/L); for the dark halo, the mass density distribution of 

' an axisymmetric isothermal sphere with a core was chosen. Hydrodynamic gas 

I ' simulations were performed for each galaxy for a sequence of five different mass 

^ • fractions of the stellar disk and for a wide range of spiral pattern speeds. These 

^ ■ two parameters mainly determine the modelled gas distribution and kinematics. 

• The agreement between the simulated and observed gas kinematics permitted us 

• to conclude that the galaxies with the highest rotation velocities tend to possess 
' very massive stellar disks which dominate the gas dynamics within the optical 

■ ' radius. In less massive galaxies, with a maximal rotation velocity of < 200kms~^, 

the mass of the dark halo at least equals the stellar mass within 2-3 disk scale 
lengths. The maximal disk stellar mass-to-light ratio in the i^'-band was found 
to lie at about M/L^ ~ 0.6. Furthermore, the gas dynamic simulations provide 
a powerful tool to accurately determine the dominant spiral pattern speed for 
galaxies, independent of a specific density wave theory. It was found that the 
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location of the corotation resonance falls into a narrow range of around three 
exponential disk scale lengths for all galaxies from the sample. The corotation 
resonance encloses the strong part of the stellar spiral in all cases. Based on 
the experience gained from this project, the use of a color-correction to account 
for local stellar population differences is strongly encouraged when properties of 
galactic disks are studied that rely on their stellar mass distributions. 

Subject headings: galaxies: kinematics and dynamics — galaxies: structure — 
galaxies: halos — galaxies: spiral 

1. INTRODUCTION 

Much evidence has accumulated in recent years that the stellar disks of high surface 
brightness spiral galaxies dominate the mass budget of the inner regions. Most of the studies 
arguing for a maximal disk scenario, are based on the detailed analysis of high resolution 
rotation curve measurements (Blais-Ouellette et al. 1999; Palunas & Williams 2000; Ratnam 
& Salucci 2000; Salucci 2001). These studies, however, generally derive the rotational support 
of the stellar disk from an axisymmctric disk mass model and allow no consideration of non- 
circular rotation components. More sophisticated modelling strategies have been applied to 
spiral galaxies with variable significance (Erickson, Gottcsman, & Hunter 1999; Pignatelli et 
al. 2001). For the most part those studies also support heavy disks for high surface brightness 
galaxies, but they also find candidates for which lower disk mass fractions are more likely. 
The most convincing direct arguments for the maximal disk scenario come from the studies 
of strong bars in spiral galaxies. These features induce a strong dynamic trace in the velocity 
field and provide a good laboratory for estimating the stellar mass component. Based on fiuid 
dynamic modelling, a maximal disk solution is found by Englmaier & Gerhard (1999) for 
the Milky Way and by Weiner, Sellwood, & Williams (2001) for NGC 4123. Furthermore, 
theoretical considerations hint for the requirement of a non-massive halo contribution in 
the central regions of strongly barred galaxies with a very cold disk, because otherwise 
dynamical friction would slow down the bar very quickly, leading to its destruction Debattista 
& Sellwood (1998, 2000). On the other hand, for initially hotter disks, this effect is largely 
reduced Athanassoula (2002) and references therein. 

However, there is also evidence that even high surface brightness spiral galaxies might 
be dominated by the dark matter mass component in their central disk regions as suggested 
by high resolution simulations of cosmological dark matter halos (Moore et al. 1999a; Ghigna 
et al. 2000; Fukushige & Makino 2001, e.g.). Bottema (1997) inferred from stellar velocity 
dispersions in spiral galaxy disks that a more massive halo component is needed to explain 
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the findings. There are recent studies that argue for hghter disk models by making use of 
previously hardly exploited relations. Courteau & Rix (1999) applied a statistical TuUy- 
Fisher relation analysis to a large sample of galaxies trying to relate the maximal rotation 
velocity of a galaxy to its disk size. Mailer, Simard, & Guhathakurta (2000), and more 
recently Trott & Webster (2002) used the geometry of gravitationally lensed systems to 
disentangle the effects of the stellar disk and the halo masses. These groups found that the 
dark halo also dynamically plays an important role in the galaxies' inner regions. 

In the present study, wc quantify the relative mass fractions of the stellar disk and the 
inner dark halo by analyzing the non-axisymmctric force component in the gravitational force 
field that arises in response to the spiral structure of the observable stellar mass distribution. 
We address this issue by exploring the influence of the spiral arms on the gas velocity fleld 
that builds up in response to the total gravitational potential of the galaxy. We model the 
gravitational potential from a mass map of the stellar distribution and a dark halo model. 
With the use of hydrodynamical gas simulations, we predict the gas velocity fields for different 
empirical mass models of the sample galaxies, assembled from a variety of stellar disk and 
dark matter mass components. The two parameters that dominate the resulting forces in 
the disk are the stellar disk mass contribution and the pattern speed of the stellar spiral, Qp. 
Thus, in this study we explore the effects of those parameters and compare the simulated 
gas morphology and velocity wiggles with our observations, enabling the identification of 
the best fitting scenario for each individual galaxy. Comparable studies have been applied 
to barred spiral galaxies taking advantage of the strong dynamical trace of the bar in the 
velocity field (Weiner, Sellwood, & WiUiams 2001). 

First results from our study have been presented already in Kranz, Slyz, & Rix (2001), 
(hereafter, KSROl) for NGC 4254 where it was proven that this strategy yields valuable 
information which allows us to constrain disk and halo fractions also for spiral galaxies with- 
out a strong bar. In the present paper we only briefiy discuss the modelling and comparison 
strategies involved (Section 2) but instead focus on the conclusions of the analysis. More 
details about the modelling issues can be found in KSROl, Slyz, Kranz, & Rix (2003) and 
Kranz (2002). 

2. OBSERVING AND MODELLING THE SAMPLE GALAXIES 

For the studied galaxies we collected near-infrared (NIR) photometry as the basis to 
study the luminous mass distribution and gas emission line (Ha) spectroscopy to acquire 
kinematic information for the systems. The NIR images were taken during two observing 
runs in May 1999 and March 2000 with the Omega Prime camera at the Calar Alto 3.5 
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m's prime focus (Bizenberger et al. 1998). It provides a field of view of 6^76 x 6'76 with a 
resolution of 0'.'3961 per pixel. We used the K'filter which has a central wavelength of 2.12 \mi. 
We obtained the gas kinematics from longslit spectroscopy measurements of the Ha emission. 
With the setup we used, the TWIN spectrograph that is located at the casscgrain focus of 
the Calar Alto 3.5 m telescope achieved a spectral resolution of 0.54 A per detector pixel, 
which translates to 24.8 km s^^ LOS- velocity resolution per pixel, allowing us to determine 
LOS- velocities with ~ 7kms"^ precision. We sampled the velocity field of the galaxy along 
14 or 16 slit position angles. The sht of the TWIN spectrograph measures 4'x I'.'S on the sky. 
The spectra were obtained during three observing runs in June 1999, May and December 
2000. 

The modelling procedure involved four separate steps. First, a stellar mass map had to 
be derived for each galaxy from the NIR images, as the basis to calculate the gravitational 
disk potential. Then, a parameterized dark matter component had to be added to the stellar 
mass component, constrained to render the observed rotation curve of the galaxy. In a third 
step wc carried out two-dimensional hydrodynamical gas simulations in these gravitational 
potentials, modelling the gas density distribution and the gas velocity field. Finally, these 
models needed to be compared to the observations to draw conclusions about the galaxies' 
structure. This detailed analysis was performed for the four galaxies NGC 3810, NGC 3893, 
NGC 5676 and NGC 6643, following KSROl. We included NGC 4254 in the final analysis, 
which was pubhshed in KSROl. 



2.1. Constructing the stellar mass maps 

For the first step, converting the stellar light map into a stellar surface mass map of 
unknown stellar M/L, we used color-corrected NIR photometry, applying a correction recipe 
provided by Bell Sz de Jong (2001). Bell & de Jong showed, based on population synthesis 
models and the assumption of a universal initial mass function, that a galaxy's color is 
tightly correlated to its mass-to-hght ratio. This apphes to a galaxy's global color as well as 
to local colors, and is approximately true, regardless of whether the color changes are due 
to dust, age, or metallicity differences. The color-M/L relation can be employed to improve 
the estimate of the stellar surface mass density distribution, S^, of galaxies taken from the 
measured NIR surface brightness, hk, by using color maps of the galaxy: 

2.5\ogJ:^ ^-liK + aK{V-K)+CK (1) 

for the example of {V — K) color maps. The optical images that were required to determine 
the galaxy colors were taken from the literature (see Table 1). Although this correction is 
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physically motivated, it relies on empirical constants a and determined from numerical 
studies. Following Bell & de Jong (2001, Table 1), we used uk — 0.314 and values for ck, 
depending on the galaxies' physical distances. 

When applied to our sample galaxies, this color correction (Eqn. 1) had two major 
consequences. First, the radial disk scale length shortened from the K-band image to the 

color corrected stellar mass image. Second, the detailed two dimensional distribution changed 
by attenuating or enhancing different azimuthal light features in the disk. Figure 1 illustrates 
how the color correction affects the radial light distribution of NGC 6643. The steepening of 
the radial profile is due to the overall color gradients within the disk. Interestingly, the color 
correction also diminishes the deviation of the stellar profile from an exponential disk, which 
in the stellar light appear in part caused by different stellar populations. Despite the color 
correction, the derived mass map still may bear some systematic uncertainties for NGC 5676 
and NGC 6643, because these galaxies are classified as starburst galaxies. The correction 
(Eqn. 1) need not be a good approximation in this regime. 

The color correction not only applies to population differences, but also works for mod- 
erate dust absorption (ropt < 1)- Dust makes the stellar radiation dimmer and redder. To 
first order, this is accounted for by Eqn. 1, with similar constants. Eventually, even in 
the 7^-band, which is usually considered as quite a good estimate of the stellar mass dis- 
tribution, the color correction changes the radial scale length of a galaxy by about 10 %. 
Consequently, this causes the inferred radial force components in the disk to change their 
magnitude, affecting the shape of the rotation curve. The use of more than only two colors 
for the correction, could further improve the quality of the mass density distribution that 
can be derived. Besides this color correction, several other minor procedures (e.g., spatial 
filtering) have been applied to the K-hand images. They were described in more detail in 
KSROl and Kranz (2002). 



2.2. Gravitational potentials and hydrodynamic models 

In order to derive total gravitational potentials, $tot, for the galaxies, a maximal stellar 
M/L was chosen for the disk, which corresponds to a maximal stellar potential contribution. 

For five different values of the stellar M/L, parameterized by a disk fraction fj, the total 
gravitational potential of the galaxy was assembled by adding the potential of a dark halo 



^The constant c combines both empirical parameters and ones that are specific to the particular galaxy 
that is studied: c = c/ + ^q, where c/ is the empirical constant and Hq is the surface brightness of 1 Mq at 
the galaxy's physical distance. 
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component, constrained to match the observed rotation curve of the galaxy. For the dark 
halo, the mass density distribution of an axisymmetric isothermal sphere with a core was 
chosen. Its gravitational potential, $haio) was also calculated in the two-dimensional plane 
of each galaxy's disk: 

^tot{R\U) = fd MR) + $halo(ii|fd) (2) 

with ffj = 0.2,0.45,0.6,0.85 and 1. Finally, a fixed pattern speed was chosen to represent 
the dominant spiral feature, leading to centrifugal and coriolis forces in the simulations. 

The hydrodynamic simulations were performed using the Bhatnagar-Gross-Krook (BGK) 
hydro-code, of Adrianne Slyz. This is an Eulerian, grid based hydrodynamics code which 
is derived from gas kinetic theory. Its concept and performance is well documented in Xu 
(1998), Slyz (1998), Slyz & Prcndcrgast (1999), Slyz et al. (2002). The gas density distribu- 
tions and gas velocity fields resulting from these simulations were compared to the observed 
galaxy morphologies and to the measured H II kinematics (see Kranz 2002 for modelling 
and comparison details). The evaluation of the data match let us probe both the relative 
luminous-to-dark matter fractions for the sample galaxies and the pattern speed of the spiral 
structure. 



3. RESULTS FROM THE ANALYSIS 

For the discussion of the results we also include the findings for the spiral galaxy NGC 
4254 which was studied with the same approach and has been presented in KSROl. In hght 
of this, our sample consists of five high surface brightness spiral galaxies for which we discuss 
the modelling results in a common context. 

3.1. Disk dynamics 

The secular evolution of spiral structures in disk galaxies is not yet fully understood. The 
debate is about whether spirals arise from quasi-stationary, self-sustained density waves, or if 
they are transient features, triggered by internal or external processes, that fade away unless 
continuously replenished (Athanassoula 1984, for a review). The present study is not very 
dependent on the long-term evolution of the spiral. The gravitational potentials used for the 
modelling are calculated from the morphological "snapshot" that the measurements provide 
and the potentials were not allowed to evolve during the simulation process. To achieve a 
steady state flow in the simulations we needed to adopt a spatially and temporally constant 
pattern speed for the dominant spiral structure. The good qualitative match of observed 
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and simulated morphology in our sample suggests that in the end this approximation might 
not be far off. Furthermore, it must be noted that the modelling only requires that the spiral 
morphology is stable over about one dynamical period, the time required for the gas to adjust 
to moderate changes in the pattern speed. Aside from some subtle evidence that weak central 
bars of our sample galaxies reveal pattern speeds that are decoupled from the rest of the 
disk, the rigid spiral pattern assumption seems to be consistent for most of the disks studied 
here. To avoid confusion with possible effects of these weak central bars, we exclude in such 
cases the very inner parts of the disk from the comparison with the simulations and thus 
restrict ourselves to the galaxies' outer disks. 



3.1.1. Location of the corotation resonance 

If a fixed pattern is a reasonable approximation, the determination of resonance locations 
within galactic disks is well defined. It has already been proven that such an approach yields 
valuable information about the disk dynamics of several galaxies, e.g., M 51 Garcia-Burillo, 
Combes, & Gerin (1993), M 100 Garcia-Burillo, Sempere, & Combes (1994), NGC 7479 
(Sempere, Combes, & Casoh 1995) and M 94 and NGC 3310 (Mulder & Combes 1996). 
However, since these studies were mostly dealing with single systems, the results were not 
put into a more general context. 

KSROl showed that it is possible with our modelling approach to locate the corotation 
resonances with a fairly high precision, by comparing the simulated gas density distributions 
for different spiral pattern speeds to the observed disk morphology. The decision for the 
best pattern speed, or equivalently the corotation radius, -RcRi was based on a visual match 
in the morphologies, which accounted for the phase pitch angle and the extent of the spiral 
arms, as well as for bifurcations. This method appears to be of a qualitative nature, and it 
might be argued that the pitch angle of the response depends on the velocity dispersion of 
the gas, in the sense that cold material will form a more tightly wound spiral than the stars. 
However, comparing the simulated gas response to images of the galaxies taken in optical 
wavelengths rather than the NIR images, it can be seen that the most dense, shocked regions 
in the simulations come to lie in fact inside the stellar spirals and tend to coincide with dust 
lanes, in cases there are some. Furthermore, this approach was supported by the fact that 
the simulated velocity field for the selected Rcr also matches best with the observed rotation 
curves and less good for other pattern speed assumtions. In this way it is also possible to 
estimate an uncertainty range. Details of the Rqr estimates are given in KSROl and Kranz 
(2002). Table 2 compiles the results for the five analyzed galaxies and Figure 2 illustrates 
the location of the corotation (circle) with its uncertainty range. Figure 3 compares the 
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estimated corotation radii to the disk scale length. Taking the X-band exponential disk 
scale length, the mean ratio of the corotation radius and the disk scale length shows little 
variation among the objects: 2.71 ± 0.43, increasing to 3.04 ± 0.49 for the color corrected 
exponential disk scale lengths. Already from this rather small sample it is possible to outline 
general trends for spiral galaxies. 

First, one can examine the spiral morphologies and how they reflect the dynamical state 
of the disk. In the modal theory of spiral structure (Bertin ct al. 1989a,b) the corotation 
resonance is typically located in the outer parts of the optical disk, where the gas-to-star 
density ratio is larger. The stellar spiral should exhibit a basic continuity across the corota- 
tion resonance and a reversal of the properties of the spiral tracers, such as the locations of 
H Il-regions (Bertin 1993). As is apparent from Figure 2, these predictions are only partly 
satisfied for the galaxies in our sample. On the other hand, some of the galaxies reveal multi- 
arm or flocculent spiral properties that cannot be explained by linear, quasi-stationary modal 
theories. Non-linear mechanisms or transient amplification processes might play important 
roles and affect the location of the corotation resonance, too. The results from non-linear 
orbital modelling of spiral structures (Patsis, Contopoulos, & Grosb0l 1991) indicate that 
the strong, symmetric, logarithmic stellar spirals already end at the inner 4:1 resonance, 
located well inside the corotation resonance. In their approach, the end of the strong spiral 
pattern is found at a radial range of 2 - 3 exponential disk scale lengths, which is compatible 
with a corotation radius of > 3 i?exp- As seen from Figure 3, the corotation radii from our 
galaxy sample tend to fall in that range. 

Figure 2 compares these corotation radii to the image morphologies. Clearly, the dom- 
inant two-arm pattern in these galaxies lies mostly within the corotation resonance. This 
finding does not disagree with the conclusions of Patsis, Contopoulos, & Grosb0l (1991) and 
Patsis & Kaufmann (1999), derived from non-hnear orbital modelling and N-body simula- 
tions respectively. However, the linear modal theory of spiral structure also predicts that 
the corotation radius falls into the radial range, which was found for the sample galaxies 
discussed here. Bertin (1993) states Rqr ~ 3 — 4i?cxp- Thus, the locations of the corotation 
resonances in the sample galaxies arc mostly in agreement with the predictions from linear as 
well as non-linear theories. The good agreement between the measured disk properties and 
the theoretical expectations from spiral density wave theories can be regarded as support of 
these theories. 

Determining the location of the corotation resonance from hydrodynamic modelling is 
a rather elaborate and time consuming process. However, looking at Figure 2 it becomes 
clear that there is no obvious common feature in the galactic disks that would enable an 
easy identification of the corotation resonance otherwise. It is not always the case that at 
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corotation the star forming rate is largely reduced. Furthermore, there is no characteristic 
feature in the amplitudes of the low Fourier components at the radius of the corotation 
resonance. This lack of any reliable corotation tracer complicates the application of density 
wave theories to real galaxies. Ultimately, the appearance of any real galaxy in the universe 
is the product of a combination of many processes. In some cases, the linear modal density 
wave theory can be applied to describe the relevant issues, but in many galaxies non-linear 
effects in the density wave dynamics and other processes seem to play a role. Nevertheless, 
assuming a fixed spiral pattern rotation speed seems, to first order, to be a viable assumption 
for the simulations of gas flows in fixed galaxy potentials. In light of this, it is a very promising 
conclusion that the hydrodynamic gas simulations provide a powerful tool to learn about 
disk dynamics and to determine the location of resonances in spiral galaxies. 



3.1.2. Quantitative match of the 2D gas kinematics 

In addition to the good agreement of the simulated gas density distribution with the disk 
morphology, the modelled gas velocity fields provide a good overall match to the measured 
kinematics too. This comparison can be made along the 14 or 16 position angles along which 
we have long-slit Ha velocities. Figure 4 gives an example of how good the simulations 
reproduce individual features in the observed kinematics. Displayed are the rotation curves 
along three slit positions across the disks of the galaxies NGC 3893, NGC 5676, and NGC 
6643. Overlaid are the simulation results from each a maximum disk and a minimum disk 
model with the best fitting model indicated by the thick line. In all cases, the global structure 
of the velocity field could be reproduced with satisfying accuracy by the model velocity field. 
This is a non-trivial result which speaks highly for the quality of the hydrodynamics code. 

In order to compare the amplitude of individual wiggles and rotation curve features to 
the observations, the code was required not only to reproduce global gas dynamics but also 
to properly model the small scale features. In general, this is also quite well achieved by the 
simulations, as it can be read from Figure 4. Slyz, Kranz, & Rix (2003) discuss the details 
of how the positions and amplitudes of features in the simulations depend on our modelling 
assumptions. However, within the galactic disks, there are processes other than local gravity 
fluctuations which create additional gas-dynamic effects. The gas dynamic traces of expand- 
ing supernova gas shells or outflows from young star forming regions are superimposed on 
the larger scale response to the gravitational potential wells and may eventually cover up 
the features that are relevant for this project. In order to minimize the influence of these 
non-gravitational effects, we used the hydrodynamic simulations - for which the only driving 
force is gravity - as a reference to discard a fraction of 30 - 40% from the observed data 
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points that deviate largely from all simulated cases as non-gravitationally induced wiggles. 
This procedure is explained in more detail in §4.3.2.1. of Kranz (2002). It should be added 
that with this treatment only local parts of the rotation curves, such as individual wiggles 
or features near the center, where small bars could play a role, were excluded and that the 
global coverage of the disk is not sacrificed. Ultimately, it would be desirable to use high 
sensitivity two dimensional velocity maps to provide a better understanding of the dynamical 
processes in the disk and to completely ehminate this concern. 

Using this technique we find that the agreement between the simulated two dimensional 
velocity fields and the observed rotation curves is very good also on a wiggle-by- wiggle basis. 
As an example, this is shown for the averaged rotation curve of NGC 3893 in Figure 5. 
About 30 % of the data points have been discarded. Although the match of the large wiggle 
at about 12" radius is not perfect, the simulated velocity field provides a better general 
representation of the rotation curve features than any axisymmetric rotation curve model 
could do. 



3.2. Ratio of stellar to dark mass 

Our main objective was to find the relative stellar disk to dark halo mass fractions 
from the comparison of modelled gas flows within realistic galaxy gravitational potentials 
of differing dark and luminous matter composition to observations. The results from the 
model-to-data comparison is presented in Figure 6 for the flve galaxies NGC 3810, NGC 
3893, NGC 4254, NGC 5676 and NGC 6643 drawing on a least squares comparison of 
predicted velocity fields for differing f^ to the observed. The symbol sizes represent the 
probability - derived from the least squares analysis - that a certain mass model applies 
to the specific galaxy. The most probable values for the disk mass fraction is indicated by 
the dotted regions in the plot. Since in some cases a purely formal least squares analysis 
does not give an unambiguous indication of the mass composition of the galaxy, we have to 
reexamine the quality of the fit by eye in order to assign a certain disk mass fraction to the 
specific galaxy. In the cases of NGC 3893 and NGC 5676 this approach supported higher 
disk mass fraction models (Kranz 2002). As can be seen from Figure 6, the precision with 
which a certain scenario can be attributed to a particular galaxy is not equally good for 
all sample galaxies. For NGC 3893 the disk mass fraction could be estimated quite well, 
placing this galaxy into the realm of maximum disks. For NGC 3810 and NGC 6643 the 
least squares comparison mostly argued for a light to medium disk scenario. The analysis 
for NGC 6643 was especially difficult since this galaxy exhibits only a weak stellar spiral 
density wave, leading to a larger uncertainty range for the disk mass fraction. The analysis 
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of NGC 4254 suffered from the high abundance of strong non-gravitationally induced gas 
kinematic features in the observed velocity field, which is apparent from the high abundance 
of small scale wiggles in the obsevcd rotation curves. The main conclusion for NGC 4254 was 
to exclude the strictly maximal disk solution. For NGC 5676 the x^-analysis mostly argues 
for a heavy disk scenario, although some simulations suffered from premature terminations. 
However, as discussed in Kranz (2002), the simulations proceeded to a point where the 
spiral pattern seemed to have reached a quasi-stationary state, making us confident that the 
conclusions for this galaxy are quite reliable. 

As defined by Sackctt (1997), the designation "maximal disk" applies to a galactic disk 
when 85 ± 10% of the total rotational support of a galaxy at a radius R = 2.2 Rexp is 
contributed by the stellar disk mass component. In the extreme case of v^isk — 0.75 Vtot, this 
definition allows a relatively massive halo to almost comprise the same amount of mass as 
the disk within that radius. 

The above definition translates into the notation used here in the following way. The 
rotational support from the stellar disk is i'disk(-R2.2) = VidV-k{R2.2), where t'*(i?2.2) is the 
rotational support from the maximal stellar disk at a radius of 2.2 disk scale lengths. A 
maximal disk model also includes a halo contribution i'*(i?2.2) = fmaxi'tot(^2.2), with fmax 
typically 0.9. Accordingly, the rotational support from the disk is about 0.9 (fd = 100%), 
0.83 (fd = 85%), 0.7 (fd = 60%), 0.6 (fd = 45%), and 0.4 (fd = 20%). All of the fd = 85% 
models clearly represent maximal disk scenarios. The fd = 60 % models are at the lower 
border of the maximal disk realm, representing compositions with about an equal amount 
of dark and luminous mass within 2-3 exponential disk scale lengths. In light of this. 
Figure 6 shows that a maximal disk scenario applies to NGC 3893 and NGC 5676. For the 
other galaxies the results from the modeUing favor less massive stellar disks. NGC 3810, 
NGC 4254 and NGC 6643 possess stellar disks that appear to balance the mass of the dark 
halo inside the optical disk (fd ^ 60 %) , although for NGC 6643 an even hghter stellar mass 
component could not be excluded. 

Table 4 gives an overview of the stellar mass-to-light ratios that were derived from the 
calculations of the potentials. The values are given for a maximum disk and mostly fall in 
the range between M/L^ = 0.6 - 0.7. A smaller value for M/Lk was only found for NGC 
4254. These stellar mass-to-light ratios agree very well with those emerging from the color 
correction procedure of Bell & de Jong (2001), which are based on population synthesis 
models^. Using these mass-to-light ratios, the total mass of the stellar disk can be inferred. 



^Bell & de Jong (2001) used an initial mass function that was scaled to a maximal disk scenario, which was 
established from measured disk kinematics. In light of this, the population synthesis M/Ls also ultimately 
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The values given in Table 4 apply to a disk mass fraction which, as taken from Figure 7, is 
a reasonable assumption for the particular galaxy. These disk masses appear realistic and 
scale well with the maximal velocities in the rotation curves, i.e. attributing a high disk 
mass to NGC 5676 and a light disk to NGC 3810. The mismatch for NGC 4254 most likely 
originates from residual inaccuracies in the determination of the galaxy's distance, inclination 
and possible calibration errors due to non-photometric conditions during the data acquisition. 
The other results for NGC 4254 are only insignificantly affected by possible inaccuracies of 
these parameters. 

4. DISCUSSION 

We have presented a new, direct and independent estimate of the stellar mass fraction 
in a small sample of spiral galaxies. With respect to ratio of luminous and dark matter in 
their inner parts, the population of high surface brightness spiral galaxies does not seem 
to comprise an entirely homogeneous class of objects. Combining the maximum rotation 
velocity, presented in Table 1, with the estimates for the galaxies' most probable fractions of 
the stellar disk mass, reveals an intriguing trend: the most massive of the analyzed galaxies, 
NGC 3893 and NGC 5676, tend to also possess the most massive stellar disks which dominate 
the dynamics of the central regions. The other galaxies from the sample are less massive 
systems that exhibit a maximum rotation velocity fc<200kms~^. For those objects the 
total mass of the dark halo within the optical radius is higher and was found to, at least, 
equal the stellar mass. This trend is graphically presented in Figure 7. 

Our results can be compared to Athanassoula, Bosma, & Papaioannou (1987). These 
authors inferred disk and halo mass fractions from spiral structure constraints. The small 
points displayed in Figure 7 represent a subsample of the total sample (given in Athanassoula 
et al. 1987, Table A2), selected for late type galaxies with a bulge mass Mbuige < 20% Mjisk 
and a maximum rotation velocity of more than lOOkms"^. For comparison the stated disk- 
to-halo mass ratio was converted into our scheme of quantifying the disk mass fraction f^. 

It can be clearly seen that the trend that emerged from our small dataset is supported 
by the results from Athanassoula, Bosma, & Papaioannou (1987). In this sample there is no 
late type spiral galaxy with a maximum rotation velocity of over 200kms~^ that seems to 
be in agreement with a substantially submaximal stellar disk. For less massive systems this 
is not the case anymore. Apparently galaxies with a wide range of dark halo and disk mass 
decompositions can exist and maintain spiral structure. 



rely on disk dynamics. 
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In light of this, the results from Courteau & Rix (1999) might be explained that in a sta- 
tistical sample of spiral galaxies the heavy and clearly maximum disk galaxies could account 
for only a relatively small fraction and most of the "lighter" high surface brightness galaxies 
exhibit already a considerable dark mass fraction. Indeed, in the Courteau & Rix sample, 
only roughly a third of the galaxies exhibit maximal rotation velocities Vc > 200kms~^. 

However, beyond the discussion of the disk mass fraction in spiral galaxies, at present 
there is growing evidence that galaxy size dark halos are not cuspy in their centers, as claimed 
by standard cold dark matter (CDM) models. (Moore 1994; Flores & Primack 1994; Moore 
ct al. 1999b; Salucci & Burkcrt 2000; Borriello & Salucci 2001; de Blok ct al. 2001; van den 
Bosch & Swaters 2001). This circumstance was first found for low surface brightness and 
dwarf galaxies, that supposedly are comprised of a massive dark halo. Most of the authors 
argued that this also applies to high surface brightness galaxies. For the present analysis, 
the functional form of the dark halo profile has only a second order effect on the results as 
long as the models describe well the overall shape of the observed rotation curves. Thus, we 
decided not to distinguish between different dark halo profiles. However, the results from 
this study, which mostly confirm the strong dynamic infiuence of the stellar disk in massive 
high surface brightness spiral galaxies within the optical radius, provide very little room for 
a strongly cusped halo core. 



5. SUMMARY 

With this project it has been demonstrated that an analysis similar to the one applied 
to barred galaxies (Englmaier & Gerhard 1999; Weiner, Sellwood, & Williams 2001) can 
also be applied to non-barred spiral galaxies. This application is on one hand more difficult 
and problematic, since the spiral arms provide a much less prominent dynamical feature in 
the velocity field than a strong bar does. On the other hand, for pure disk galaxies the 
vertical scale height is smaller in the central regions than in strong bars that often show 
peanut shapes when viewed edge-on. This fact reduces systematic uncertainties that arise 
from motions along the z-axis which are superimposed onto the circular motions that are 
studied. The main points that could be learned from this study can be summarized in the 
following way: 

• The color correction, described in Section 2.1, allows us to derive the stellar mass 
density from any multi color imaging data. As a correction to the K-band photometry 
Equation 1 has proven to work very well. We encourage its use for all studies that 
depend on the mass distribution in galaxy disks. 
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• The BGK hydrodynamics code, apphed to the problem (Slyz, Kranz, & Rix 2003), is 
a powerful tool and it can be used successfully to simulate the interstellar gas in the 
disks of galaxies. The parameters that dominate the morphology and the kinematics of 
the gas simulations are the non-axisymmetric stellar mass component and the pattern 
speed of the spiral structure. 

• Adopting a globally fixed pattern speed for the rotation of the spiral seems to be a viable 
approximation for the simulation process. Moreover, from analyzing the morphology of 
the simulated gas density distribution it could be found for the sample of five galaxies 
that the corotation radius Rcr = 2.71 ±0.43i?cxpi increasing to 3.04 ± 0.49-Rexp for the 
color corrected exponential scale lengths representing the stellar disk mass. In all five 
cases the strong stellar spiral ends inside the corotation resonance. 

• The detailed comparison of the velocity fields is challenging, but it is possible to put 
constraints on the relative mass fractions of the inner disk. The results from this study 
indicate that high surface brightness galaxies possess maximal disks if their maximal 
rotation velocity Vc > 200kms~^. If the maximal rotation velocity is less, the galaxies 
appear to have submaximal disks. 

We would like to thank the Calar Alto staff for their support in obtaining the data for 
this study. AS acknowledges the support of a fellowship from the UK Astrophysical Fluids 
Facility (UKAFF). 
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Radius in arcsec 

Fig. 1. — The effects of tlie color correction, displayed for tlie example of NGC 6643. Dis- 
played are tlie B-band, K-band and color corrected profiles. Due to the outward blueing of 
the disk, the color corrected profile becomes steeper than the other two profiles. The smooth 
lines are fits with double exponential profiles. 
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Fig. 2. — Location of the corotation resonance, -Rcr; including its uncertainty range, over- 
laid onto deprojected, spiral enhanced i^-band galaxy images. The -Rcr's were derived from 
comparing the spiral arm morphology with the gas morphology in hydrodynamical simula- 
tions of differing pattern speeds. The values for i?cR are listed in Table 2. Note that the 
corotation resonance appears to lie just outside the strong spiral. 
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Fig. 3. — Location of the corotation radius in terms of the exponential disk scale length for 
the sample galaxies. The asterisks mark the values for the i^-band disk, while the boxes 
refer to the color corrected scale lengths. The error bars come from uncertainties in the 
corotation radius estimations. The horizontal lines represent mean values; the lower one for 
the i^'-band disk, the upper one for the color corrected disk. 
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Fig. 4. — Examples of the comparison between simulation results and observed kinematics 
(data points) for the galaxies NGC 3893, NGC 5676, and NGC 6643. Displayed are three 
position angles (PA) for each of the galaxies with velocity fields from "maximal disk" and 
"minimal disk" simulations overlaid. The best fitting model is indicated by the thick fine. 
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Fig. 5. — Effect of tlie two dimensional simulation on the rotation curve (averaged around 
the major axis), displayed for NGC 3893. The two dimensional fit (continuous line) provides 
a better representation of the inner rotation curve than the axisymmetric fit (dash-dotted 
line) does. 
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N3810 N3893 N4254 N5676 N6643 

Fig. 6. — Results from the least squares analysis. Larger symbols represent better x^-values. 
The symbol sizes are scaled to give equal sizes for the best fitting points. Regions of best 
agreement between simulations and observations are marked by the dotted regions. Note 
that these shaded areas also account for a visual inspection of the fit quality. 
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Fig. 7. — The dependence of the disk mass of galaxies on their total mass, as inferred from 
the maximal rotation velocity. The disk mass fraction f^ is defined in Equation (2). In 
addition to the results from this study (filled dots) results from the study of Athanassoula, 
Bosma, & Papaioannou (1987) (small dots) are included in the plot. The galaxies can be 
separated into two classes. The most massive galaxies tend to possess disks with high mass 
fractions, f(j. Less massive systems may exhibit a considerable dark mass fraction. 
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Table 1. Properties of the sample galaxies. 



Galaxy 


Btot 


K'tot 


PA 


i 


max Vj-ot 


dist. 


optical 


reference 




[mag] 


[mag] 


[1 


n 


[kms"-*-] 


[Mpc] 


image 




NGC 3810 


11.4 


8.02 


22 


46 


155 


13.5 


B 


Frei et al. 1996 


NGC 3893 


11.23 


7.94 


166 


42 


220 


17 


B 


Tully et al. 1996 


NGC 4254 


10.2 


6.85 


67.5 


41.2 


153 


20 


g 


Frei et al. 1996 


NGC 5676 


11.7 


8.03 


46.5 


63.2 


245 


33 


V 


Heraudeau & Simien 1996 


NGC 6643 


11.8 


8.18 


39 


57.5 


185 


23 


B 


Heraudeau & Simien 1996 
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Table 2. Corotation radii and exponential scale lengths for the galaxies. 



Galaxy Rcr Re^{K') -Rexp(i^'corr) 

[kpc] [kpc] [kpc] 



NGC 3810 3.15 ± 0.5 1.07 0.913 

NGC 3893 5.5 ± 0.5 1.80 1.74 

NGC 4254 7.5 ± 1.1 3.54 3.06 

NGC 5676 11.0 3.59 3.12 

NGC 6643 6.5 tli 2.72 2.48 



Note. — All values are calculated for the distance 
given in Table 1. -Rcxp(-^'corr) specifies the color cor- 
rected exponential scale lengths. 
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Table 3. x^/A'"-values from the least squares analysis. 



Galaxy 


fd = 20 


fd = 45 


fd = 60 


fd = 85 


fd = 100 


NGC 3810 


0.800 


0.550 


0.555 


0.787 


1.069 


NGC 3893 


1.775 


1.555 


1.480 


1.897 


1.451 


NGC 4254 


1.146 


0.959 


1.007 


1.206 


2.303 


NGC 5676 


0.706 


1.109 


1.817 


0.908 


1.463 


NGC 6643 


1.400 


1.349 


1.651 


2.284 


2.183 
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Table 4. Maximum disk stellar i^-band mass-to-light ratios and the associated stellar disk 

mass. 



Galaxy 




(M/L)eolor 


-^disk 1 max 


Mdisk 1 fd 


NGC 3810 


0.63 


0.55 


1.48 X 10^° 


0.89 X 10^° 


NGC 3893 


0.56 


0.54 


2.32 X 10^° 


2.32 X 10^° 


NGC 4254 


0.23 


0.74 


3.53 X 10^° 


2.12 X 10^° 


NGC 5676 


0.67 


0.83 


9.55 X 10^° 


8.12 X 10^° 


NGC 6643 


0.71 


0.65 


4.23 X 10^° 


2.54 X 10^° 



Note. — Given are the maximal disk K-hand M/Ls. (M/Li^)jnax 
corresponds to a stellar mass fraction fd = 100; (M/L)coior were de- 
rived using the relations by Bell & de Jong (2001) for the respective 
overall galaxy colors. The disk mass Mjisk | fj has been scaled for the 
galaxies' most probable disk mass fraction fj. All values are given 
in solar units. 



